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the iron region are studied. Allowed and first- forbidden contributions to the cross 
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sections are calculated within the quasiparticle random phase approximation, ex- 
tended to finite temperatures within the Thermo-Field-Dynamics formalism. The 
GTo strength distribution at finite temperatures is calculated for the sample nucleus 
54 Fe. The neutral-current neutrino-nucleus inelastic cross section is calculated for 
relevant temperatures during the supernova core collapse. The thermal population 
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of the excited states significantly enhances the cross section at low neutrino energies. 
In agreement with studies using a large scale shell-model approach the enhancement 
is mainly due to neutrino up-scattering at finite temperatures. 

PACS: 26.50. +x; 23.40.-s; 21.60.Jz; 24.10Pa 

I. INTRODUCTION 

Neutrinos play a decisive role in core-collapse supernova explosions since they carry most 
of the gravitational binding energy released. The transport of neutrinos through the hot and 
dense stellar environment is believed to ultimately be responsible for a successful explosion, 
although the details are not fully understood yet. The present paper addresses the role of 
thermal effects in the inelastic neutrino-nucleus scattering in the iron core during infall and 
shortly after bounce. 

At the end of 1980th it was pointed out by W. C. Haxton that inelastic neutrino-nucleus 
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scattering (INNS) mediated by the neutral-current can be of importance comparable with 
the other processes of neutrino down-scattering The INNS contributes to the neutrino 
opacities and thermalization during the collapse phase, the revival of the stalled shock 
wave in the delayed explosion mechanism, and to explosive nucleosynthesis. The estimates 
by Haxton were based on nuclei in their respective ground states, i.e. for a "cold" nuclei. 
Subsequently, it was realized that the INNS occurs in hot stellar environment (T > 0.8 MeV) 



e changes of the INNS 
2j and then in ^| on the 



and, due to the thermal population of nuclear excited states, sizea : 
cross section are to be expected. The effect was firstly analyzed in 

basis of large-scale shell-model (LSSM) calculations. In Refs. (si, it was found that the 
INNS cross section noticeably increases at T ^ and for neutrino energies E v < 10 MeV, 
especially for neutrino scattering off even-even nuclides. 

However, in the subsequent core-collapse supernova simulations jlj] including several 
dozens of nuclides, it was demonstrated that the inclusion of the INNS process does not 
have a large effect on the collapse dynamics and the shock wave propagation. But it signif- 
icantly modifies the spectrum of neutrinos generated in the v e burst. 

Here, we apply an alternative approach for treating the thermal effects for INNS cross 
sections. In essence, our approach is based on the thermal quasiparticle random phase 
approximation (TQRPA). We apply it in the context of Thermo-Field-Dynamics (TFD), 
which enables a transparent treatment of thermal excitation and de-excitation processes 
and offers the possibility for systematic improvements. This approach has recently been 
used in studies of the electron capture on hot iron and germanium nuclei under stellar 

n 

conditions 6 . 



II. FORMALISM 

A. Fundamentals of the Thermo-Field-Dynamics 

Thermo-Field-Dynamics [t 9 ] is a real-time formalism for treating thermal effects in 
quantum field theory and non-relativistic many-body theories. The standard TFD formalism 
treats a many-body system in thermal equilibrium with a heat bath and a particle reservoir 
in the grand canonical ensemble. The thermal average of a given operator A is calculated as 
the expectation value in a specially constructed, temperature-dependent state |0(T)) which 
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is termed the thermal vacuum. This expectation value is equal to the usual grand canonical 
average of A. In this sense, the thermal vacuum describes the hot system in the thermal 
equilibrium. 

To construct the state |0(T)), a formal doubling of the system degrees of freedom is 
introduced. In TFD, a tilde conjugate operator A - acting in the independent Hilbert space 
- is associated with A, in accordance with properly formulated tilde conjugation rules [7H9]. 
For a system governed by the Hamiltonian H at T = 0, the whole Hilbert space at T 7^ is 
spanned by the direct product of the eigenstates of H (H\n) = E n \n)) and those of the tilde 
Hamiltonian H having the same eigenvalues (H\n) = E n \n)). The important point is that, 
in the doubled Hilbert space, the time-translation operator is not the initial Hamiltonian H, 
but instead the thermal Hamiltonian % = H — H . This implies that the excitations of the 
thermal system are obtained by the diagonalization of H. 

The thermal vacuum is the zero-energy eigenstate of the thermal Hamiltonian T-L and 
satisfies the thermal state condition 7 

A\0(T)) = ae w/2T l t |0(T)), (1) 

where a = 1 for bosonic A and a — i for fermionic A. 

As it follows from the definition of H, each of its eigenstates with positive energy has the 
counterpart - the tilde-conjugate eigenstate - with negative but the same absolute energy 
value. This allows to treat excitation- and de-excitation processes at finite temperatures. 

Obviously, in most practical cases one cannot diagonalize % exactly. Usually, one resorts 
to certain approximations such as Hartree-Fock-Bogoliubov mean field theory (HFB) and 



the Random-Phase Approximation (RPA) (see e.g. 10|]). In what follows the TFD studies 



for neutrino induced charge-neutral excitations in hot nuclei are based in part on the results 



of 
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1^ (see also [6|). 



B. Charge-neutral excitations in hot nuclei 

In what follows we employ the Hamiltonian of the Quasiparticle-Phonon Model (QPM) 
Hqpu 13] which consists of proton and neutron mean fields H sp , the BCS pairing interactions 
impair and isoscalar and isovector separable particle-hole interactions. Since the inelastic 
neutrino-nucleus scattering involves nuclear J n excitations of both natural (7r = (— 1) J ) and 
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unnatural (ir = (— 1) J+1 ) parities both the separable multipole H M h and spin-multipole H$ M 
interactions are included in the particle-hole channel 

Hqpm = H sp + impair + + H$ M . (2) 

The four terms of Hqpm read 

H sp = E E ( E 3 ~ ^r) a ]m a jm ' 
T=p,n jm 

#pair — E ^ T E a }m a jm a y^T a j'm'y 

r=p,n jm 



j'm' 



H m =~E E (4 A) +P^ A) )Mi(r)M v (pr) . 

A/i rp=±l 

< = ~E E (4 LA) +P«i LA) )^(r)5 LAM (pr), 
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LA/i rp=±l 

Here, we use standard notations of the QPM. Namely, a] m and a jrn are the creation and 
annihilation operators of particle with quantum numbers jm = n, l,j, m and energy Ej\ jm 
stands for the time reversed single-particle states; the index r is isotopic one and changing 
the sign of r means changing n ■<->■ p ; the parameter G T is the constant of pairing interaction; 
A T is the chemical potential; the parameters (k^) denote the strength parameters of 
the isoscalar (isovector) multipole (a = A is a multipole index) and spin-multipole (a = LA 
is a spin-multipole index) forces. The multipole M£ (t) and spin-multipole S^ Xfl (r) single- 
particle operators read as 



n m 2 



M U T ) = E T 0'i m il^^A(r)^b2m 2 )a] imi a j2 

jlTOl 

h m 2 

«Uw = E T 0'i m ii iL ^( r )[ y ^Ki^)< mi s 2m2 , (3) 



h m 2 



where 



[ y H£ = E( LMlm i A ^) y ^» CT - > 

M,m 

and the notation ^ r implies a summation over neutron (r = n) or proton (r = p) single- 
particle states only. The excitations of natural parity are generated by the multipole and 
spin-multipole L = A interactions, while the spin-multipole interactions with L — A ± 1 are 
responsible for the states of unnatural parity. 
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To determine the thermal behavior of a nucleus governed by the Hamiltonian (J2J) we 
should diagonalize the thermal Hamiltonian "Hqpm = -£^qpm — Hqpm and find the corre- 
sponding thermal vacuum state. This will be done in two steps. 

In a first step, the sum of single-particle and pairing terms "Hbcs — %s P + "Hpair is di- 
agonalized. To this end two subsequent unitary transformations are made. The first is 
the usual Bogoliubov u,v transformation from the original particle operators aj m , a- m to 
the quasiparticle ones aj m , a jm . The same transformation is applied to the tilde operators 
a] m , a.- m , thus producing the tilde quasiparticle operators 8?j m , S jni . The second, unitary 
thermal Bogoliubov transformation mixes the original and tilde degrees of freedom 

fijm = X j a ]m~Wj a jm (4) 
0jm = ''/»./„. • '!/;",». (tf + V 2 = !)• 

The operators 0j m ,/3j m ,(3j m , and (3j m are called thermal quasiparticle operators. 

The coefficients found by diagonalizing "Hbcs an d demanding that the 

vacuum of thermal quasiparticles is the thermal vacuum in the BCS approximation, i.e., it 
obeys the thermal state condition (CQ). As a result one obtains the following equations for 

Uj, Vj and Xj, yf 

V] = ^(l-^±y\u 3 = {l-v^l\ (5) 
y^fl + exp^)]" 172 ^^^-^ 2 ) 172 , (6) 



where Ej = a/ {Ej — A T ) 2 + A 2 . The coefficients y 2 determine the average number of ther- 
mally excited Bogoliubov quasiparticles in the BCS thermal vacuum 

(0(T); qp|aj m a im |0(T); qp) = y) (7) 

and, thus, coincide with the thermal occupation factors of the Fermi-Dirac statistics. 

The pairing gap A r and the chemical potential A T are the solutions to the finite- 
temperature BCS equations 

^( T ) = ^ L E r ^ + 1 )( 1 - 2 ^>^' 

N T = ^(2j + l)(^ 2 x 2 + u)y% (8) 



6 



where N T is the number of neutrons or protons in a nucleus. 
At this stage, the thermal BCS Hamiltonian Hbcs is diagonal 



and corresponds to a system of non-interacting thermal quasiparticles. The vacuum for 
thermal quasiparticles |0(T);qp) is the thermal vacuum in the BCS approximation. The 
states (3j m \0(T); qp) have positive excitation energies whereas the corresponding tilde-states 
(3j m \0(T); qp) have negative energies. Since the thermal vacuum contains a certain number 
of Bogoliubov quasiparticles, excited states can be built on |0(T);qp) by either adding or 
removing a Bogoliubov quasiparticle. The first process corresponds to the creation of a 
non-tilde thermal quasiparticle with positive energy, whereas the second process creates a 
tilde quasiparticle with negative energy. 

At the second step of the approximate diagonalization of "Hqpm, long-range correla- 
tions due to the particle-hole interaction are taken into account within the thermal QRPA 
(TQRPA). Within the TFD formalism the terms "H^ 1 and "Hg^ are written in terms of the 
thermal quasiparticle operators determined above. Then, Hqpm is approximately diagonal- 
ized within a basis of thermal phonon operators 
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+ C-J/v/y; + + *i?&\Mh\}\> ( 9 ) 



where [ denotes the coupling of single-particle angular momenta ji, J2 to a total angular 
momentum A. Now the thermal equilibrium state is treated as the vacuum |0(T);ph) for 
the thermal phonon annihilation operators. 

The thermal phonon operators are considered as bosonic ones which imposes certain 
constraint on the phonon amplitudes. To find the amplitudes and energies of the thermal 
phonons, the variational principle is used, i.e., we find the minimum of the average value of 
thermal Hamiltonian with respect to the one-phonon states Q^|0(T); ph) or Q1_^.|0(T); ph) 
under the aforementioned constraint. 

After variation one obtains a system of linear equations for the amplitudes V^Vo Ww,; Vh^i 
etc. as well as for the energies (details can be found in ref. These constitute the 

equations for the thermal quasiparticle random phase approximation. In contrast to the zero 
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temperature case, the negative solutions of the secular equation have a physical meaning. 
They correspond to the tilde thermal one-phonon states and arise from f3* terms in the 
thermal phonon operator. As it was noted above, creation of a tilde thermal quasiparticle 
corresponds to the annihilation of a thermally excited Bogoliubov quasiparticle. Conse- 
quently, excitations of negative-energy thermal phonons correspond to transitions from 
thermally excited nuclear states. 

After diagonalization in terms of thermal phonon operators the TQRPA part of the 
takes the form 

?W = ^(Ql^Q^ - Ql^Q^). (10) 

Xfii 

To fix properly the thermal vacuum state |0(T);ph) corresponding to TRPA we once again 
turn to the thermal state condition (Cp) and derive the final expressions for the amplitudes 
of the thermal phonon operator (]9]). 

Once the structure of thermal phonons is determined, one can determine the transition 
probabilities from the thermal vacuum to thermal one-phonon states. They are given by the 
squared reduced matrix elements of the corresponding transition operator 7a^ 

^H^ilTAiio^ph)! 2 , 

SA^KgAiiiTAiiocr^ph)! 2 . (ii) 

Thus, the probability to excite the hot nucleus is given by $Ai, while $Ai is the probability 
to de-excite it. 



C. Cross section of inelastic neutrino-nucleus scattering 

Considering neutrino-nucleus inelastic scattering in stellar environments we assume that 
a nucleus is in thermal equilibrium with a heat bath and particle reservoir or, in TFD terms, 
in the thermal (phonon) vacuum state. An inelastic collision of a hot nucleus with neutrinos 
leads to transitions from the thermal vacuum to thermal one-phonon states. 

In the derivation of the relevant cross section at finite temperature we follow the formalism 



by Walecka-Donnelly 14], 115], which describes in a unified way electromagnetic and weak 
semileptonic processes by taking advantage of the multipole decomposition of the relevant 
hadronic current density operator. In the case of neutral-current neutrino-nucleus scattering, 
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the differential cross section for a transition from an initial nuclear state (i) to a final state 
(/) can be written as a sum over all allowed multipolar it ies J n 



da. 



2G 2 {E v -u if ) 2 cos 



2 oo oo 

^{£^+£<4 



7T 



2J,; 



J=0 



J=l 



where 



and 



(12) 



(13) 



e 



i(^ii^r s ii^)r+i(^ii^ii^)r_ 

2Re(J f \\J™^J l )(J f \\jf\\J i y 



e 



e 



- tan 2\/"2^ + tan2 2 



(14) 



Here G is the electroweak coupling constant, O is the scattering angle, E v is the incoming 
neutrino energy, Uif is the transition energy from the initial nuclear state (i) to the final 
state (/), and q li = (tUif,q) (q = \o\ = y ojff + AE V {E V — oj^) sin 2 ® j is the four- momentum 
transfer. The operators Mj, Lj, Jf, and J™ ag are the multipole operators for the charge, 
longitudinal, and the transverse electric and magnetic parts of the four-current, respectively. 
Following 14J they can be written in terms of one-body operators in the nuclear many-body 
Hilbert space. 

The cross section involves the reduced matrix elements of these operators between the 
initial and final nuclear states. Within the present approach, the initial nuclear state is 
the thermal phonon vacuum (TV) and the final states are the thermal one-phonon states. 
Therefore, at T ^ all the reduced matrix elements in Eqs. ( TT3lfT4|) are calculated in 
accordance with Eqs. (fTTj) . The total cross section is obtained from the differential cross 
sections by summing over all possible one-phonon states of different multipolarity and by 
numerical integration over scattering angles 

1 da^v^j 



/e{Ai} ' 



dVt 



-dcosQ. 



(15) 



Up to moderate energies [E v ~ 15 — 20 MeV), the inelastic neutrino- nucleus scattering 
is dominated by the neutral-channel Gamow- Teller transitions J n = 1 + . Moreover, in the 
q — t- limit, the full operator exciting 1 + states is reduced to the following Gamow- Teller 
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operator: 



GT = (— ) at (16) 



where {qa/9v) = —1-2599 16] is the ratio of the axial and vector weak coupling constants, 
a is the spin operator and to is the zero-component of the isospin operator in spherical 
coordinates. 

To circumvent computational limitations in the LSSM calculations Ufl the total INNS 
cross section a(E u ) was split into two parts - a down-scattering part a^{E v ) and the up- 
scattering part a u (E u ). The term <7d(E u ) includes transitions where the scattered neutrino 
loses energy whereas the term a u (E v ) includes those transitions where the neutrino gains 
energy from a hot nucleus. Assuming the validity of the Brink hypothesis for the GT 
resonance, the down-scattering term was transformed to a sum over only those final excited 
nuclear states which are coupled by a direct GTo transition with the nuclear ground state. 
As a result, <Jd(E u ) appeared to be independent of T. 

In our case, the part aa(E u ) corresponds to transitions from |0(T); ph) to \Qxi) states with 
positive energies whereas the a u (E u ) term is the sum of transitions |0(T); ph) — > \Qm) where 
the tilde-states have negative energies. In the latter transitions a neutrino gains energy due 
to nuclear de-excitation. 

Thus within the present approach the GT (J 7 " = 1 + ) contribution to the cross section 
reads 



a(E u ) = a d (E u ) + cr u (E u ) = — V(£, - Uji) 2 $ji + — Y t (E v + u M ) 2 <$> Ju (17) 



The probabilities $jj and $jj are given in (TIT]) with T = GT . Since u>ji, $j« and $jj are 
functions of T, both terms a d and a u depend on temperature. 

Whereas the GT component determines the neutrino- nucleus cross section at low E u , 
higher multipole contributions become increasingly important at higher neutrino energies. 
Moreover, at higher neutrino energies Eq. (fl6l) for GTq is not valid and the 1 + transition 



operator will depend on transfer momentum q. According to Refs. 
reduces the cross section. 
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181 ] the (/-dependence 
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III. CALCULATIONS FOR THE HOT NUCLEUS 54 FE 

Numerical calculations have been performed for 54 Fe. The single-particle wave functions 
and energies were calculated in a spherically symmetric Woods-Saxon potential. The con- 
stants of the pairing interaction were determined to reproduce experimental pairing energies 
in the BCS approximation. All parameters are the same as in our previous calculations 
jf], for electron capture rates on the same nucleus at T ^ 0. 

The radial dependence of the residual multipole and spin-multipole forces is chosen in the 
form R\(r) = dU(r)/dr where U(r) is the central part of the single-particle Woods-Saxon 
potential. Thus, R\(r) as well as the parameters Kq X 1 and do not depend on A. The 
isovector parameters and k[ lx ^ are fitted to the experimental position of the El 20 1 



and Ml [2l| resonances in 54 Fe. According to the estimates in Refs. 22|, |23|, the isoscalar 
spin-multipole interaction is very weak in comparison with the isovector one. Following 23], 
we take Kq 1 ^ / = 0.1. 

First, we have performed TQRPA calculations of the GTo strength distribution in 54 Fe. 
As in the LSSM calculations [4|, the GTo operator (|T6l) have been scaled by a quenching 
factor 0.74. In Fig. [TJ we display the GT strength distributions for the ground state 
(T = 0) of 54 Fe and at three stellar temperature values, occurring at different collapse stages: 
T = 0.86 MeV corresponds to the condition in the core of a presupernova model for a 15M 
star; T = 1.29 MeV and T = 1.72 MeV relate approximately to the neutrino trapping 
and neutrino thermalization stages, respectively. All results are plotted as a function of 
the energy transfer to 54 Fe. For charge-neutral reactions this energy is equal to a thermal 
phonon energy Uji. 

At T = 0, the transition strength is concentrated mostly in one-phonon 1 + state forming 
the GT resonance near u ~ 10 MeV. The main contribution to the phonon structure comes 
from the proton and neutron single-particle transitions I/7/2 —> 1/5/2- With temperature 
increase the fraction of low-energy transitions in the GT strength distribution increases. 
The physical reason is the weakening and subsequent collapse of pairing correlations (at 
T pa 0.8 MeV) and appearance of low-energy particle-particle and hole-hole transitions due 
to thermal smearing of neutron and proton Fermi surfaces. Moreover, at finite temperature 
the "negative energy" transitions to tilde one-phonon states appear. As a result, the GT 
energy centroid is shifted down by 1.1 MeV at T = 1.72 MeV. This indicates a violation of 
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the Brink hypothesis within the present approach. 

The contribution of 1 + transitions to the INNS cross section is shown in Fig. [2]^a) for 
different temperatures. The calculations have been performed with the exact g-dependent 1 + 
multipole transition operator [l^j]. As in the LSSM calculations 3|, the cross section <j(E v ) 
at T = is equal to zero when E v is less than the energy of the lowest 1 + state in 54 Fe. 
Within the QRPA, the lowest 1 + state in 54 Fe has an excitation energy of w(l + ) ~ 7.5 MeV 
(see Fig. CQ. The GT transitions at T ^ do not show such a gap due to thermally 
unblocked low- and negative-energy transitions. As a consequence, there is no a threshold 
energy for neutrinos at finite temperatures and the INNS cross section appears to be quite 
sensitive to T at neutrino energies E u < 10 MeV. As it follows from the present calculations 
as well as from the LSSM study [3J], thermal effects can increase the low energy cross section 
by up to two orders of magnitude when the temperature rises from 0.86 MeV to 1.72 MeV. 
Finite temperature effects are unimportant for E v > 15 MeV where excitation of the GT 
resonance becomes possible and dominates the cross section. These features were pointed 
in as well. 

To check the influence of finite momentum transfer on the INNS cross section we also 
have performed calculations with the GT transition operator (fl6|) . A comparison of 1 + 
and GTo cross sections is shown in Fig. EJ^b) for T = 0.86 MeV. The g-dependence becomes 
important at E v > 30 MeV. At E v = 35 MeV the INNS cross section calculated with the 
g-dependent 1 + operator is by 20% less than that calculated with the GT operator ( fl6|) . 
At E v = 50 MeV the difference is by about factor of 2. The effect does not change with 
temperature. 

The contribution of first-forbidden transitions 0~, 1~, and 2~ to the INNS cross section 
were also calculated within the TQRPA, taking into account the g-dependence as given 
in [14J]. The results are presented in Fig. |3j As it can be seen, a temperature increase 
enhances the cross sections at low and moderate E v . The main reason is thermally unblocked 
low-energy first-forbidden transitions. According to our calculations 2~ transitions dominate 
the total contribution of first-forbidden transitions to the cross section at low neutrino 
energies, while at higher energies the total contribution is mainly determined by the 1~ 
transitions. 

In Fig. HJ the INNS cross sections at different temperatures are shown as a sum of 
1 + ,0~,1~, and 2 _ contributions (we omit the contribution of the + multipole because 
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it is negligible). At low E v the cross sections are almost completely dominated by the 
GT transitions. The part of the cross sections arising from the first-forbidden transitions 
becomes increasingly important at larger E v . We find that for E v = 30 MeV up to 20% 
of the cross section is due to first-forbidden transitions. For E u = 40 MeV allowed and 
forbidden transitions contribute about equally, while at E v = 50 MeV the contribution of 
first-forbidden transitions is nearly twice as large as that of 1 + transitions. 

In the LSSM calculations, the temperature-related enhancement of a(E u ) was only due 
to the neutrino up-scattering. In our approach both the up-scattering and down-scattering 
parts of a(E u ) are temperature dependent. To analyze the relative importance of these two 
types of scattering processes we display them separately as the functions of E v for different 
values of T in Fig. El 

A weak T-dependence of cx^ is seen at low neutrino energies E v < 12 MeV. At higher 
energies practically does not depend on T. As the function of E v the down-scattering 
cross section sharply increases at low neutrino energies and then grows more slowly. Instead, 
o u is quite sensitive to temperature but its dependence on E„ is obviously smoother than 
that of <Jd (at least at E v < 15 MeV). The absolute values of cr^ and o u are of the same 
order of magnitude only at quite low neutrino energies E u < 4 — 10 MeV. 

Thus the conclusion is that the T-dependence of the INNS cross section at low neutrino 
energies is mainly due to up-scattering process whereas at neutrino energies E v > 15 MeV 
when the thermal effects are much less important the INNS cross section is determined by 
the neutrino down-scattering. 

The above conclusions agree well with the results of the LSSM studies for even-even nuclei 
{3, 4|. Furthermore, our results for 0^ confirm the applicability of approximations based on 
the Brink hypothesis, which has been used in calculations of in the LSSM. 



IV. CONCLUSIONS 

We have performed studies of the temperature dependence of the cross section for inelastic 
neutrino-nucleus scattering off the hot nucleus 54 Fe. Thermal effects were treated within the 
thermal quasiparticle random phase approximation in the context of the TFD formalism. 
These studies are relevant for supernova simulations. 



In contrast to the large-scale shell-model studies 



4j we do not assume the Brink hy- 
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pothesis when treating the down-scattering component of the cross section a(E u ). Moreover, 
we take into account thermal effects not only for the allowed 1 + transitions but also for the 
first-forbidden transitions CF, 1~, and 2 _ . For all multipole contributions we have performed 
the calculations with momentum dependent multipole operators. 

Despite these differences between the two approaches, our calculations have revealed the 

nn 

same thermal effects as were found in [3J, |4j: A temperature increase leads to a considerable 
enhance of the INNS cross section for neutrino energies lower than the energy of the GTo 
resonance. This enhancement is mainly due to neutrino up-scattering at finite temperature. 
The calculated cross sections for 54 Fe are very close to those given in (4]. Thus, the results 
of our study show that the present approach provides a valuable tool for the evaluation of 
the inelastic neutrino-nucleus cross sections under stellar conditions. The approach can be 
easily adopted to calculate the INNS cross sections as a function of scattering angle. 
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Figure captions 

Fig. 1. GT strength distributions in the 54 Fe nucleus at different temperatures T as a 
function of the energy of transition u. 

Fig. 2. a Contribution of 1 + transitions to the cross section of neutrino inelastic 
scattering off 54 Fe calculated with the g-dependent 1 + excitation operator as a function 
of neutrino energy E v at different stellar temperatures T; b - A comparison of the cross 
sections of neutrino inelastic scattering off 54 Fe calculated with the g-dependent 1 + excitation 
operator (solid line) and the GT excitation operator (Tl6|) (dashed line) at T = 0.86 MeV. 

Fig. 3. Contributions of different first-forbidden transitions to the neutrino-nucleus 
inelastic scattering cross sections for 54 Fe at different temperatures T: a - the contribution 
of the 0~-transitions; b - the contribution of the 1 "-transitions; c - the contribution of 
2~-transitions; d - the summed contribution of the all first-forbidden transitions. 

Fig. 4. The neutrino-nucleus inelastic scattering cross sections as the sum of allowed 
and first-forbidden contributions for 54 Fe at different temperatures T. 

Fig. 5. The down-scattering ad(E u ) (a) and the up-scattering a u (E u ) (b) parts of the 
neutrino-nucleus inelastic scattering cross section for 54 Fe at different T. 
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